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I*  Introduetion 

Th«  algoritha  of  this  papor  ik  baaad  on  tha  dual  linaar  prograwdnf 
algorithm  ([1],  (6])and  tha  raoant^  davalopad  boond  aaoalation  aathod  for 

Tha  t^rld-dual  algorithm  la  diirldad 


into  two  atagaa.  In  tha  firat  ataga  a  linaar  prograa  ia  aolrad  to  /laid 
an  optimal  aolutlon  in  i^aetiocial-Taluad  rariablea.  Tha  final  tablaau 
giran  by  tha  dual  linaar  programaing  mathod  ia  than  tranafornad  and  axpandad 
into  a  naw  tableau  of  a  madil/  apecifiable  canonioal  form*  In  tha  aaoond 
ataga  of  tho  algorithm  a  variant  of  tha  bound  aaealatioo  mathod  ia  appliad 
to  tha  naw  tablaau— -aod  to  tha  tablaaua  auccaaaival/  darivad  tharaaftar«' 
until  ona  or  mora  of  a  diatinguiahed  aat  of  oolumna  (and  a  oorraaponding 
aat  of  rowa)  attains  a  pradataminad  configuration.  At  thia  point  tha 
indioatad  rowa  and  columns  ara  diacardad,  navar  to  ba  raoovarad,  and  tha 
mathod  oontinuaa  recuraivaly  with  tha  bound  aacalation  algorithm 
until  tha  problam  ia  solved^  _ ^ 


laooihl 


It  may  ba  obsarvad  that  tha  aaoortd  ataga  of  tha  hybrid-dual  algorithm 
worka  almost  axaotly  in  revarsa  of  Oomory's  original  intagar  programming 
mathod  [U]-  With  tha  hybrid-dual  msthod,  onca  tha  fractional  linaar 
programaing  tablaau  is  tranafornad  and  axpandad,  no  naw  oonstrainta  or 
variablaa  ara  addad  tharaaftarj  inataad  cartain  variablma  and  conatrainta 
ara  dalatad  as  tha  aolutlon  procasa  avolvaa..  In  oontraat,  tha  original 
Oomosy  msthod  bagina  with  tha  final  linaar  programming  tablaau  and 
progressivaly  incraaaaa  ita  dimenaion,  though  tha  tablaau  may  ba  kapt  from 
axeaading  a  eartain  aisa  by  waading  out  naw  variablat  whan  tbay  amarga  from 
tha  nonbaaio  sat. 


\ 
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1  Introduction 

The  algorithm  of  this  paper  ie  based  on  the  dual  linear  pro^’raniraing 
algorithm  ([1],  (6])and  the  recently  developed  bound  ercalation  method  for 
solTing  integer  linear  programs  (3j  The  hybrid-dual  algorithm  is  divided 
into  two  stages  In  the  first  stage  a  linear  program  is  solved  to  yicild 

k 

an  optimal  solution  in  fractional  valued  variables  The  final  tableau 
given  by  the  dual  linear  prograimtiing  method  is  then  transfoimed  and  expaided 
into  a  new  tableau  of  a  readily  specifiable  canonicel  form  In  tne  sec o  id 
stage  of  the  algorithm  a  variant  of  the  bound  escalatio;'.  inethod  is  appli*  d 
to  the  new  tableau-  and  to  the  tableaus  successively  derived  thereafter- 
until  one  or  more  of  a  distinguished  set  of  columns  (and  a  corresponding 
set  of  rows)  attains  a  predetermined  configuration  At  this  point  the 
indicated  rows  and  columns  are  discarded,  never  to  be  recovered,  and  the 
method  continues  recuisively  with  the  bound  escalation  algoritlun 
until  the  problem  is  solved 

It  itwy  be  observed  that  the  second  stage  ol  the  hybrid-duel  algorituri 
works  almost  exactly  in  reverse  of  Gomory  s  original  integer  progranmrn  ' 
method  II4]  With  the  hybiid-dual  method,  once  the  fraction.il  linear 
progremraing  tableau  is  transfor-iied  and  expanded,  no  new  constraints  or 
variables  are  added  thereafter}  instead  certain  variables  and  constraint r 
are  deleted  as  the  solution  process  evolves  In  contrast,  the  original 
Gomory  method  begina  with  the  final  linear  programming  tableau  and 
progressively  increases  its  dimension,  thougii  the  tableau  may  be  kept  frim 
exceeding  a  certain  size  b/  weeding  out  new  variables  when  they  emerge  from 


the  nonoasic  set 


Two  additional  featurea  of  &he  hybrid-dual  algorithm  deaerve  nentiohr 
First,  it  la  poaaibla  to  work,  with  a  aonewhat  analler  tableau  than  tho 
standard  bv  following  one  of  two  modified  prooedurea  inTolving  r'atrictwd 
rulea  of  choice.  Second,  the  hybrid -dual  algorithm  uaually  eonvergea 
before  all  of  the  problem  conatrainta  are  completely  aatiafied  When 
the  tableau  reaches  a  prescribed  form,  one  portion  of  the  bottom  row  may 
be  added  to  another  to  yield  the  optimal  aolution  imaediateLy.  In  this  <4ay 
it  is  frequently  poasibLe  to  reduce  the  number  of  steps  o'-herwise  required 
to  solve  the  problem. 

The  ner.t  section  Intioduces  the  notational  framework  to  be  used  in 

this  papar,  along  with  a  brief  discussion  of  the  principal  characteriatiis 

of  the  dual  linear  programming  algorithm  and  the  bound  escalation  method  ^ 

The  h/br'd-dual  "IgoritluD  is  outlined  in  Se  .tion  ill.  Theorems  concernl-ig 

the  properties  of  the  algorithm  are  stated  tr.id  prored,  inrludirg  tbe 

fundaxaental  reeuH  that  convergence  to  an  optimal  aolution  is  asaured  in 

a  finite  number  of  steps  for  any  bounded  problem  with  a  nonempty  solution 

set  In  Sevtlon  IV  an  oxaiqple  problem  la  solved  in  two  different  ways  to 

applied  to  the 

demonstrate  t)ie  work.'nga  of  the  hybrid-dual  *.Lgorlthm  and  c<nputatlonal 
standard 

fsatuT'^s  tableau  and  also  to  a  coi^ieased  tableau  In  conclusion  aosM 
com.-ui-a  .ional  features  of  the  mthod  are  discussed. 

1 1 .  Notation  and  C<^tituent  Algorithms 

Ihe  integer  programming  problem  mty  be  written 
Minijaire  wb  b 

c 

•ubject  to  wA  >  c, 

w  >0,  w  integer 

1,  It  is  not  strictly  necessi-ry  to  use  the  dual  linear  programming  algorithm 
as  opposed  to  the  priaial  algorithm  I/s  have  however  selected  the  former  in 
order  to  facilitate  subsequent  exposition- 


-s- 


where  b  ts  an  id*!  column  vector,  c  is  a  l>in  row  vector,  Is  a 

scalar,  A  is  an  m»*n  matrix,  and  w  is  a  l^m  row  vector  whose  com^ionent-s 
we  wis^i  to  find  in  order  to  satisf/  tne  inequalities  and  optimize  tlie 
mialAiMtion  criterion,  fte  aesuaa  Uiat  the  components  of  A,  b,  ai'd  c 
are  all  integer.  Without  further  loss  of  generality  we  also  specify  the'' 
the  augmented  matrix  (  b  A)  is  lexicographically  negative  by  row  (aee 
f2]  and  16]).  The  initial  tableau  for  the  dual  linear  progranuuinc  algorithm 
may  then  be  represented  as  follows. 


(1) 


-Vi 

A 

man 

*^nxn 

b 

o 

c, 

l^^n 

1 

The  first  column  of  the  tableau  representst  the  objective  function  ti- 
be  optimised,  and  each  subsequent  column  identified  a  problem  constraint 
obtained  by  rewriting  the  inequality  wA  >  c  in  UiO  form  wA  wl  •<  c, 

t’ 

wiiere  w*(w  .,...,w  )  ia  »  Ixn  row  vector  of  nonnegativa  "elack'' 

variables..  lYie  dual  linaar  programming  method  recuraively  trsiisforrta  the 
given  tableau  into  a  new  one  by  a  aequence  of  pivo..  reduction  operations 
At  any  step,  the  pivot  element  is  determined  by  selecting  a  positive 
entry  in  the  current  c  vector,  then  detarisining  t.e  unique  poaitive  entry 
in  the  associated  column  which  will  leave  the  up^er  purtioni  of  the  tableau 
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f.icalijr  attnr  carrying  out  th*  pivot  op*«itli»ii.^  Tiit 

proceiiu  i3  coa;  ietro  wtcn  thvr*  iire  nf>  positive  ontrlea  in  the 

'  tranaforsM-d)  r  vector  to  e^iooao  fiMuu  4  ponltive  component  at  tha  base 
ni  a  cx.X^.ii  Otiose  •utries  are  all  ric.-ative  or  zero  sigtiols  that  the 
‘  CdSihie  aoiutloi.  3f  t  of  the  pnihl<  a  it-  emwty, 

Jan  <lLkii.ation  in  the  lual  lutaar  roKraas.ing  al>rorltJn 
;jr(  -er.e!*  the  Idaritlt^  of  the  urip^inal  varlablee,  aifi  each  UibleAu  iofinait 
a  rro  :e.:  ex>fCtly  e<)ul/aleiit  to  U.c  ongliial.  ll^ncc  (1)  :.ay  be  used  to 
r^preuei.t  na  arbit,'*arj'  tableau  o'Aair.el  with  tae  dual  al.'oritha  after 
rt  rt'C  /alues  of  thi.  vurlabiee  at  the  nptital  solution  are 

th  1  obt  isied  by  netting  Uie  vbriablee  as.^oelated  Mfltii  the  final  A  iGAtrlc 
e  )uai  to  7.>-ro. 

fri'  *  tl.e  a  nn  notation,  the  li.ltial  tableau  for  the  bound  escalation 
T  t.’io-j  cw-:  tx-  eplcted  as  lo^i.awa. 


r  — . 

)  "^rjcl  j  ‘.MU 

^  ) 

1 

1 

mx.u 

i  «>  i  c, 

1  o  1  Ixn 

1  1 

IXK 

here,  as  befose,  the  first  c.iluun  recr<.aei.ts  th«:  objective  fuTiCtlon 
.3^..  the  e  .rcescive  ooluans  Identify  the  problem  conatra Lnts.  Tltls  tLne, 
►  ever,  t.if  canitraiats  art  ieft  as  .nequalJties  an  -  the  nor.r*ej^llv1t* 
ron  -Jt-nr  e  h.-'  au...  ar  ZOO  (ui.t>Iicltly  in  t!<e  Last  m  columns, 

I  I  .ot  .-lij'^it  it  always  -1,  ooWi-nea  by  oAiitiplylr.>t  tier  pivot 

Cl  lu’-u.  b.-  t  e  apfmrriate  roriStai.t,  The  pivot  operation  thaii  'aikeo  all 
At..er  enfif  s  iXi  tSc  row  as  the  pi/ot  clei:tert  equal  to  stro  by 

j  ii-'  I.  rectPSrti",  ;.iul*  1,  le  of  the  pivot  coluruj  to  each  of  tne  otner 

•'aiu  u.  i..  Si.fi't  ts lon.j 
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Because  the  bound  uactlation  mottiod  works  sooswtiat  dlf fsrsntljr  than 
uhs  linssr  progrsu'iij^g  algorithm^  It  is  eonvsnisnt  to  rsprssent  ths 
i^snarsl  form  of  the  tsblssu  by  igxwilng  ths  partition  bstwosn  ths  A  ostrix 
siki  the  iasntity  ostrix,  ylslding  s  tsblssu  of  ths  for* 


“^axl 

b 

o 

^lx(»H») 

T 


Ths  solution  to  ths  Intsgsr  progps».-dn*  proble*  is  obUinsd  with 
the  bound  ososistion  asthod  when  ail  elaosiits  of  ®jj((*+n)  **"*  *>***^P®®^^^^** 
The  optinsLL  w  Tsctor  then  appaars  as  the  negative  of  the  vector  originally 
lying  beneath  ths  1  oatrix. 

To  appl/  ths  bound  escalation  aethod  it  is  desired  to  create  a 
tableau  idiich  upon  suitable  indexing  B»y  he  partitioned  as  follows. 


(4) 


”^mxl 

P*P 

e 

®(B-p)xp 

b 

0 

(1, 

bep 

* 

Here  D  is  a  pXp  square  oatrix  with  positive  entries  along  the 
aain  ilagonal  and  nonpositive  entries  everywhere  else.  Q  Is  a  aatrix 

conposed  entirely  of  nonpositive  entries,  ano  d  la  a  p-eomponent  row 
vector  at  least  one  of  whoee  entries  is  positive.  We  ere  unconcerned  with 

the  portions  of  the  tableau  marked  with  a  atar. 


The  bound  esoslation  aetnod  operetes  on  such  e  structure  to  obtain 


a  vector  d»  satisfying  the  following,  two  conJltlors;  (1)  d  -  d»R  ie 
i.oripoa it i ve  ixiteger,  and  (ll)  each  of  the  f.rat  p  coi^xxieiits  of  the 


optinal  solution  vector  is  at  least  as  large  as  the  corresponding 

elcaaiit  of  Once  d*  is  aatemined,  the  botton  row  of  the  tableau 

is  ctiangad  bgr  subtracting  d*P  froa  it,  where  P  is  the  ovitrix  consistii^ 

of  the  first  p  rows  of  the  tableau.  Pjr  incorporating  this  procedure 

into  a  recursive  process  for  creatinK  the  subnatrices  U,  Q,  and  d,  the 

bouni  escalation  netltod  converges  to  the  optlaml  Integer  solution  in  s 

2 

finite  number  of  atepe. 

It  will  suffice  here  to  eoikSider  onl>  the  case  in  which  D,  C,  ard 
d  conatltute  a  single  coluosi  of  the  tableau.  In  this  instance,  d  and  0 
each  are  composed  of  a  single  positive  element,  d^  and  ratpcctivelj , 

and  d*  eonaists  of  the  single  elenteni  ^ 

li  the  vector  least  one  positive  componci^t*  it  i» 

easy  to  ersata  a  column  of  the  desired  fora  bj'  the  following  rules. 

First,  select  a  posit Iva  entry  from  the  bottom  row  of  the  tablaauc  All 
but  one  of  the  remaiiLing  positive  entries  in  its  column  nay  than  ba 
elljlnaled  by  consecutively  selecting  any  pair  of  them  and  subtracting  a 
positive  integer  multiple  of  the  row  in  shieb  one  appaara  from  the  row 
containing  tho  other.  The  only  retrlction  on  this  proceas  is  that  tha 
upper  portion  of  the  tableau  oatst  be  kept  lexicographically  negatlwa. 

Once  a  column  of  the  daelred  Ibrs  is  obtained,  it  is  peraissible  to  add 
an  integer  nailtiple  of  the  unique  row  with  the  positive  component  to  atiy 
other  row,  provided  no  entriea  that  were  previously  nonpositlva  (ir  the 
specified  column)  thereby  become  positire.  More  generally,  row  subatrac- 
tloiis  art  permissible  under  any  circumstances  so  long  as  the  tableau  in 

1.  The  optimal  solution  vector  ie  defined  relative  to  the  problem  represent¬ 
ed  by  the  current  tableau,  and  changes  whan  the  tableau  is  transformed. 

2.  D,  Q,  and  d  need  not  appear  explicitly  in  Ute  tablcmu,  but  may  be  generated 
by  any  nonnegative  linear  combination  of  the  tableen  columns  (except  for  the 
first) . 

3.  We  will  use  to  denote  the  smalleat  integer  greater  than  or  agual  to  x. 
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utair.talnod  Icodco^aphlcally  ncigatlve,  but  row  additions  csn  be  carried 

out  only  relative  to  Uia  D,  Q«  d  coiifi^urotion  in  the  oanner  indioated* 

3 

(ixcept  Uist  tne  comronent  of  d  is  allowed  to  be  sero.  These  basic  re¬ 
sults  provide  all  that  is  neee;>sary  for  the  eectioua  to  follow. 

Tl»e  ^L^b^..d«Dual  AlRoritha. 

We  will  let  (1)  represent  the  final  tableau  obtained  by  the  dual 

f 

linear  pro.'’rai]iiilng  algoritlua  after  reitidexin^.  Let  be  the  ositrix 

t  r  “y  w 

each  of  whose  entries  Is  given  by  a^^  ■  a^,  -  L*ijJ »  ^mxn  ^ 

the  luatrlx  with  entries  a^  »  »  »'bere  a^  denotes  the  element  In 

the  ith  row  andjth  oolumr  of  A.^  Sin<llarly«  let  e^  be  the  IXn  row  vector 
whoso  .1th  entry  vector  i<dth  Jth  entry 

w  .  I’.  I 

foAa  of  the  initial  tablsau  for  the  hybrid-dual  algorithm 

is  then  given  by 


(5) 


“'’bxI 

mxu 

I 

mxm 

a'' 

cocn 

^nxl 

-I 

-  ram 

I 

nxn 

b 

0 

o 

°lxn 

Ixn 

We  eharactorize  the  hybrid  xhial  algorithm  >s  follows. 

1*he  hybrid -dual  alRoritha, 

Stage  lo  Solve  Uie  fraetioral  linear  programolr^  problem  with  ttia 
iiuil  linear  progr.Juaing  methed*  and  create  the  transformed  tableau  (5). 


1.  The  nc’tation  used  to  indicate  the  greatest  integer  Jess  than  or 

e<:)ual  t<^  X. 

2.  It  is  shown  in  Section  IV  how  a  sotaewiat  smaller  tableau  omy  be  used 
by  observing  certain  additional  rules. 

3  It  may  in  fart  be  negative,  so  long  as  d^  -r  dj^  >  0. 
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St4jte  2>  Apply  the  bound  •acAlAtion  nethod  to  the  tmnafonDed 
tableau »  subject  to  the  following  Dodiflcations . 

A,  If  at  any  stage  an  element  of  the  portion  of  the  tablran 
le  less  than  eero,  its  column  umy  be  treated  aa  the  negative  of  itself 
except  in  the  determination  of  lexicographic  ordcrir^r 

0.  If  an  element  of  the  portion  of  the  tableau  le  equax  uo 
ero,  ■='nd  if  .ill  but  ono  of  the  elainents  of  Its  ooluan  are  .ilso  zero, 
ti.cn  the  indlce,tcd  coiam  and  tho  iXM  containing  tlic  elcuant  may 

be  romo/ed  from  the  tableau. 

C..  If  an  eleuent  of  c^  is  equal  to  sero,  and  If  more  than  ore 
component  of  itu  ccluran  is  different  from  earo,  the  column  (and  some 
associated  row)  may  be  pit  in  the  fora  for  removal  by  the  method  of 
Djeoren  3  (to  follow)  „  In  this  fashion  n  rows  and  n  colunns  of  tableau 
(5)  may  be  reraoved  ii  a  finite  number  of  nteps. 

f  w 

Do  The  problcaa  is  solved  when  the  vector  sum  c  +  c  is  non*- 

I* 

positive  into.f'er  (with  deleted  entries  of  c  set  equal  to  zero),  and 

the  portion  of  the  bottom  row  ie  nonpositive.  The  components  of 

tiw)  optimal  solution  voctor.  including  the  values  of  the  slack  variables 

f  w 

aro  then  read  from  the  bottom  row  to  the  right  of  c  ,  with  c  replaced 
by  c  +  c  „ 

To  prove  the  vali  iity  of  the  algorithm  we  present  the  following 
Theorems  arxl  proofs. 

Theorem  1 .  The  integer  p:  ogi*amniing  problem  BUjaTjarize*!  by  tableau  (1) 
nay  bo  solved  with  the  bound  escalation  method  applied  to  tableau  (5)» 
subject  to  the  quaiificalLon  instruction  Ao 
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I  roof;  Eftcb  of  the  coluor^  of  the  final  dual  tableau  (1)  repreaatitr 
an  equallt7  of  the  fom 


1-1 


"  Vj  "  for  J-1,  ...  .  n. 


.<hen  thl:>  is  factored  into  an  equivalent  pair  of  inequalitlee 
we  obtain 


i 


-  Vf 


'J 


To  repreaent  the  problem  in  a  tableau  stiltable  for  appl/lng  the 
bound  escalation  method,  we  thus  have^ 


(6) 


I 

mxl 

A 

nxn 

-A 

noto 

^(m4n'x(iirfn) 

^nxl 

“^nxn 

I 

nxn 

! 

o  I 

-C  - 

.r7\ _ 

^lx(n+n) 

K’e  asaunc  first  tliat  b  >  0.  For  the  final  dual  linear  progrcuBEilng 
tableau,  each  component  of  c  is  nenpoeitive.  heuco  c  >  0.  Relative 
to  the  columns  whose  bottoei  x*ow  (.dfuponents  lie  in  -c$  it  is  therefore 
poBbible  to  carry  out  the  row  additions  and  row  subtractions  indicated 
in  the  previous  section.  Since  b  >  0.  it  Is  permissible  to  subtract 
any  positive  integer  multiple  of  a  row  in  the  partition  of  the 

tableau  from  any  rew  ir.  the-A  paitition.  We  first  select  only  thooe 


1.  Vibilc  the  bound  escalation  C'^thod  is  guaxanteed  to  converge  for  an  al 
integer  tableau,  the  proofs  in  [  3]  also  obviously  apply  when  A.  b.  and  e 
are  rational.  Also,  since  the  original  problem  was  all  integer,  the  slac' 
variables  added  for  the  dual  linsar  orogxturdLng  .Igorithu  must  be  Integer 
as  well. 
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eleaente  of  -A  which  are  pt^ltlve,  and  subtract  a  lar^e  enough  multiple 

of  the  appropriate  I  „  isatrix  row  to  make  each  of  then  lees  then  or 

nui 

equal  to  zero.  Thus  for  each  element  -a^^  of  -A  such  that  -a^^  >  0»  it 
is  sufficient  to  subtract  <  >  times  the  Jth  row  of  leaving 

[  a^j  }-  a^j  («  -  <  -*^j>  )  in  place  of  the  original  -aj^j  enti7. 

The  effect  in  other  portions  of  the  tableau  will  be  to  replace  a^^  in  the 
A  netriji  by  -  C  (for  j  <  0)»  *nd  to  replace  the  0  in  row  1 

and  colujai  m  +  J  of  the  ^(gf+n)x(m+n)  ^  f  •  '•^®n  thie 

procese  is  carried  to  completion  the  onljr  positive  elements  in  the  col- 
tmins  above  the  -c  vector  will  be  those  in  the  matx*lx. 

At  this  point,  it  'uay  be  potsibie  to  add  integer  crultiplee  of  the 
iviws  of  I  V  ^  some  of  the  rows  of  (the  new)  >A.  Since  all  elements  of 

lULn 

-A  nairt  remain  nonpoeitive  after  the  addition,  the  largest  multiple  of 
the  Jth  row  of  that  can  be  added  to  the  ith  row  of  -A  (for  ^ 
is  3  •  This  quantity  will  be  zero  for  all  elAents  of  -A  that  were 

irdtiall/  nor^egative,  so  that  tlie  multiple  may  be  defined  entirely  in 
ter^  of  the  initial  -A  aatidx  (for  The  adjustmente  throughout 

the  tableau  can  thus  be  specified  in  exactly  the  oane  way  as  for  >  Oo 
If  ■■  0  the  Indicated  aujustm'nts  apply  trivially. 

The  eftanges  in  tho  c  vector  prcsczrlbod  by  the  bound  esealatioa 
method  are  readily  determined  by  the  remarics  of  the  previous  section. 

For  each  column  J  associated  with  an  element  of  d^  ■■  -c^  and 
d^  »  1,  hence  d*  "  <  •  Subtracting  the  appropriate  multiplea  of 

the  rows  of  the  I  .  matrix  from  the  vector  yields  the  value 

nXii 

(  CjJ  -  Cj  (•  -Cj  -<  )  fo**  entry  cf  the  original  -c. 

Carrying  out  the  subtraction  for  the  entire  bottom  row  of  the  tableau 
replaces  each  c^  of  the  c  vector  with  ]  ,  and  replaces  the 


-  u  - 


(t  +  J)-th  coniponeiit  Oa  the  vector  by  [c^  1,  all  other  detents 

rer^lnlng  unchanged. 

Since  ell  cosiponents  of  tne  -c  portion  of  the  tableau  are  now  non- 
positive,  we  will  not  i:  mediately  be  concerned  with  them  in  applying  the 
bound  escalation  method.  I.o.oover,  it  ie  evicent  tnat  the  colurns  assoc¬ 
iated  with  t.'.e8e  components  vdll  always  be  the  negetive  of  tr.oee  columns 
associated  with  tne  original  c  vector.  Hence  we  may  use  the  letter  to 
surauirixe  the  fort.er,  provided  we  observe  that  each  cf  these  columns  cAy 
be  treated  as  its  negative  (positioned  suitably  to  the  right ).^  Collap¬ 
sing  the  tableau  as  indicated,  we  obtain  tableau  (5)* 

Kote:  To  maintain  strict  lcxico^_rapbiic  negativity,  if  the  final  tableau 
obtained  with  the  dual  linear  programming  algorit.ja  is  degenerate— tnat 
ie,  if  some  component  of  the  fL-al  b  vector  is  equal  to  zero— t.*.en  the  forn 
of  the  6tart.ing  tableau  for  tr.e  S'Kond  stage  cf  the  hybrld-diisl  al  jorithj; 
uAj  nave  to  oe  modifird  slightly.  Perhaps  the  siiUTilest  way  to  accomplish 
tne  necessary  changes  is  to  begin  with  tableau  (5)  at  a  slightly  more 
prit-ative  stSee.  Thut,  if  bj^  »  0,  and  if  is  the  first  nonzero  coip- 
porent  in  the  kth  row  of  A,  we  define  afj  “  ®iJ*  0  (for  i  "I,,..,!  ), 

cj  •»  Cj,  cj  "  0.  Ail  other  intriee  of  A^,  A^,  c^,  c''  are  i^'lver  at 

before.  In  practice,  'vhen  o  /  0  tne  minor  departure  from  strict  leecico- 
grapnic  negativity  whicn  ray  arise  from  using  tableau  (5)  as  originally 
defined  ie  not  likely  to  hiampcr  tne  convergence  process — and  may  in  fact, 
speed  it  up. 

rheorem  2t  hows  and  colu-imis  of  tne  tableau  deleted  according  to  instruc- 
tion  b  are  irrelevant  for  obtaining  tne  optical  integer  solution. 

1.  The  "negative”  columns  nay  in  fact  be  ignored  in  the  determination  cf 
lexicographic  orderir.g  since  we  nay  consider  tnem  appended  at  the  far  id.g!'t 
of  t'bleau  (5).  The  linear  independence  of  the  rows  of  (5)  assures  that 
nons  of  c^ese  rows  will  becoL.e  all  zero. 
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Proof;  Tach  coluKn  of  the  tableau  aseociated  trith  siur-carlaes  an  «xact 
equality  in  nonne^^ative  variablse,  hence  a  deleted  ccliiicn  represents  an 
expression  of  the  form  az  »  0.  The  variable  x  associated  with  tne 
deleted  row  r.tust  ce  equal  to  zero,  and  hence  the  row  Is  superfluous. 

Cnee  the  rw  is  reiroved,  the  remainder  of  the  coiuim  is  obviously  supez^ 
fluous  as  well. 

r 

Theorem  3.  li  an  elar.ent  of  c  is  equal  to  zero,  then  the  following 
li.ethod  will  bring  the  assiciated  .column  into  the  required  form  for  rec«valr 
enabling  ti.e  deletion  of  n  rows  and  n  colunns  of  (5)  in  a  finite  nuirbor  of 
stepsr^ 

Ly  the  procedure  outlined  in  Section  II,  put  the  Indicated  column  in 
the  fem  1),  ;i,  d.  Ferforci  any  necessary  row  additions  with  the  row  conta:  n- 
ing  (cf  D)  until,  for  each  el<  aent  q  ofQ,  +  q>  0.  If  all  element  , 
oi  are  equaJ  to  zero,  the  process  is  c(xcpleti.d.  Otherwise,  consider 
negative  of  t:,e  coliran  Just  aerivod,  and  repeat. 

Proof;  '.-.e  will  use  a  superscript  to  denote  tne  step  on  which  a  given  D, 

c,  stnicture  is  oubC.ned  in  the  above  procedure.  Wo  note  that  (the  value 

•  f  d,  or.  the  secorvl  etepjnaist  be  lorived  irom  the  components  of  <^2^  .  Sli  ce 

I  j.its  of  tiic  positive  elei.ents  of  can  be  increased  ty  the  row  subtractions 

2 

soecifitc  Ir.  Section  II  for  creating  the  D,  4,  d  form,  it  follows  that  ^ 
q*,  where  -q^  is  the  r.axliaut  coc.ponant  of  But  the  .'nethod  specifiee 

t.nat  dj  +  q'^  >  0,  neiiie  >  d^.  StaJlarly,  d^  >  d^^  +  h,  whox^  h 

is  a  fixed  positive  quantity.  Thus  the  elements  of  -Q  (Which  are  non* 
r.'^gat.ive  and  strictly  Isss  then  ;:.U8t  eventually  all  be  driven  to  sero 
If  it  is  possible  for  iroLleo  tc  be  solved  >  in  widch  caso  the  elenwnt  3  cf 
c  oan  be  forced  to  zcr*o  -  all  a  of  the  corresponuiqg  colunns  (and  sofae  :i 
may  therefore  eventually  be  el  jlnated. 

1.  If  all  nonzero  cjcpone.''ts  *>f  the  colucn  have  the  same  si  71,  then  by  the 
reasoning  in  the  proof  of  Th»  reo  2.  a"*!  rows  cf  tho  tableau  containing  th.e 
norizero  elements  nay  be  eliminated  simultaneously.  We  note  also  that  the 
optlnal  solution  may  ^e  obtained  before  the  indicated  n  rows  and  eolusns  are 

ael&tad 


Tbeoren  h  Tfie  optimal  solution  vector  to  the  problem  is  obtained  at  in 
ine  traction  D^, 

Proof;  By  the  nature  of  the  bound  eEcalatlon  method,  the  vector  to  the 

right  of  (initially  the  vector  of  (6)"^  gives  integer  values 

for  the  original  problem  variables  «4hich  will  produce  an  objective  function 

value  equal  to  the  current  b^,  and  under8<)tj.8iy  (or  oversatisfy)  each  of 

the  original  constraints  by  the  positive  (or  negative!  amount  of  the 

corresponding  cor^jonent  of  the  current  vector-  At  any  stage  we  may 

increment  or  decrement  the  values  of  the  variables  associated  with  the 

original  notrlx  of  (6)  without  changing  the  objective  function 

value «  Suppose  thst  it  is  possible  to  adjust  tte  variable^  iB  cueh 

•  'jgy  ^.hat  (a)  all  of  the  problei  constraints  are  satisfied,  and  (b)  all  c  T 

the  problon  variables  are  nonnegative  integer.  SlXiCe  b^  la  monotone  non- 

decreasing  from  one  tableau  to  tt;e  next,  it  follows  that  tlie  solution  so 

obtained  is  optimal..  But  to  satisfy  the  first  n  constraints  as  strict 

equalities  (if  no  other  varirblea  are  to  be  changed),  each  of  the  '1^^ 

variables  mst  be  net  equal  to  the  corresponding  component  of  ■  This 

f  w 

nay  be  represented  by  adding  c  to  c  «  The  conditions  of  instruction  D 
then  correspond  to  t);e  conditions  for  achieving  optimality,  and  the  theoru-a 
Is  proveda 

n 

We  may  note  parenthetically  that,  since  the  components  of  c  are 
always  integer,  the  conpenents  of  c^  st  the  optimal  solution  must  be 
integer  also  Moreover,  if  lexicographic  ordering  of  the  tableau  ia 

f 

rigidly  maintained,  tha  first  nonee.''o  ca^>onent  of  c  must  be  negative,  rjnea 
with  the  bound  escalation  matnoc  the  bottom  row  of  the  tableau  must  be 
strictly  increasing  ^  lexjcofpraphically)  each  time  it  ia  changed 


Theorem  5  The  l\ybrid-»l>ial  elgcrithm  will  converga  In  a  finite  number 
of  steps  for  any  problem  containing  a  nomoptj  integer  solution  set 
bounded  for  optimality :• 

Proof  t  The  theorem  follows  innediately  from  Theorems  1  through  It  and 
the  con'/ergonce  properties  of  the  dual  linear  programming  algorithm  and 
the  bounc  escalation  method > 

iV .  ExBByle  Problem  nd  Computational  gaqperienee  ^ 

To  illustrate  the  hybrid  dual  algorittua  we  will  now  solvs  the  folloM 
ing  problem 


Minimize 
subje't  to 


^1  ♦  iWj  ♦  Itw^ 


-Iw^  «  PWj  ♦ 


1-3  >  7 


•  ‘>'2  *  *5  i  * 

"i-  "3  ;  ° 


stage  I  (the  dual  linear  programming  algorithm)! 

(o) 


0 

-1 

? 

1 

(i)  -7/2 

1/2 

5/2 

1/2 

3 

2 

1 

-1 

0 

-1 

0 

0 

-It 

1 

2 

3 

-S/2 

-1/^ 

5/2 

7/2 

U 

-1 

0 

0 

-3/? 

1/2 

1/2 

.1/2 

0 

0 

-1 

0 

0 

0 

-1 

0 

0 

_ 0^_ 

0 

_-l_ 

0 

0 

0 

-1 

0 

1 

8 

8 

21/2 

•7/2 

9/2 

25/2 

-22/7 

U/7 

16/7 

-1/7 

0 

-1 

0 

0 

0 

0 

0 

.1 

-ld/7 

5/7 

5/7 

1/7 

0 

0 

-1 

0 

-5/7 

-1/7 

5/7 

2/7 

151/7 

as/-* 

-5/7 

-25/7 
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The  optinal  solution  road  from  the  final  tabloau  Is  thus  ■  0,  •  13/7, 

Wj  -  2i/? ,  ■  d,  v.^  •  3/7,  ••  0,  yielding  an  objective  fun''tlon  value 

of  131/7. 


Tl.e  sjjqple  form  of  the  bound  escalation  taethod  outlined  in  Section  IT. 
will  suffice  for  the  second  stagao  We  will  indicate  by  an  arriw  the  colam 
relative  to  which  the  rcM  additicns  and  subtractions  are  to  be  racda  For 
definltenoas,  all  such  operatiors  w.t.11  be  carried  out  with  the  lexico 
graphically  least  negative  row,  chosen  from  among  those  in  which  the 
indicated  column  has  a  positive  entry, ^ 

We  obtain  the  initial  tableau  for  Stage  2  oy  putting  the  coiif)onents 
of  the  w  vector  in  the  order  (Wj^,  w^,  w^,  w^,  w^.,  w^). 

Stags  ? 


(o) 


(1) 


-22/1 

V7 

2/1 

<>/7 

1 

0 

0 

0 

2 

1 

-15/7 

3/7 

1 

5/7 

l/’7 

0 

1 

0 

0 

0 

0 

•5/7 

0 

0 

_1 

— -L.. 

0 

0_ 

u 

i 

0 

,  . 

0 

c 

0 

0 

1 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

1 

0 

_ 9 

0 

0 _ 

1 _ 

_  0_ 

0 

_0 

_ 0^_ 

_  0 _ 

1 

131/7 

1/7 

a/7 

5/7 

0 

u 

0 

-2 

■♦1 

•a 

1 

•1 

-2 

i 

0 

1 

0 

-3 

3 

2 

-1 

6/7 

-3/7 

2/7 

-1/7 

0 

j 

-1 

1 

0 

0 

5/7 

V? 

3/7 

2/7 

0 

0 

1 

1 

0 

0 

0 

-1 

J 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

X 

0 

-15/7 

_ 

9/7 

-1/7 

0 

3 

_-3_ 

0 

1 

11»6/7 

-17/7 

*5/7 

-1/7 

0 

0 

-3 

1 

-1 

u 

r  This  rule  is  not  always  the  most  efficient  for  either  the  simple  or  the 
■ore  general  form  of  the  bound  escalation  aeUiod  With  a  slight  embellish- 
■ent  the  rule  provides  the  baeia  for  Gomor/'s  all  integer  algorithm  [5], 
which  can  be  showh  to  be  a  variation  of  the  bound  escalation  mathod  under 
restricted  slternstives  of  choics  (see  {5}  ) 


In  o'ottixiSs^i  the  naxt  tablMui,  tho  fourth  coluan  beeanas  all  sero  except 
in  the  second  row,  hence  the  indicated  row  and  eoltan  are  deleted. 

(ii) 


-1 

-2 

1 

'  1 

0 

-3 

3 

2 

-1 

-3 

0 

1 

0 

2 

1 

1 

0 

0 

0 

•1 

A 

0 

0 

0 

1 

0 

0 

0 

0 

'•1 

0 

0 

0 

0 

1 

0 

1 

0 

-1 

4 

0 

1 

22 

-2 

M. 

•1 

0 

*1 

-2 

0 

-1 

t 


The  problsB  ie  now  solTsd.  'Wia  add  the  final  (--S  ~1  O)  to  the  final 
e'  (O  '1  -A),  yielding  the  optieid  aolntion  Teotor  (0  '•!  *2  ^ 

Wo  vlll  not  stop,  however,  but  vdU  coitimie  with  Stafa  2  of  the  alforitla 
to  illustrate  the  tx^ese  of  eliainatinf  row*  and  ooliana  from  the  tableam- 


(lU) 


■  1 

-1 

1 

0 

~3  1 

1 

(M 

~3 

1 

0 

2 

-1  1 

0 

0 

c  ! 

1 

1 

0 

0 

0  0 

1 

0 

1 

-1 ; 

1 

0 

-1 

4  -1 

0 

1 

1 

22  { 

0 

•*1 

2  -2 

■“1  'A 

' 


(iT) 


H  the  optimal  integer  eolation  the  vector  (0  *1  -2  *i2  •A) 


in  the  final  tableam  oorreeponde  to  the  nogative  of  the  reordered  vector 
(v^,  w^,  w^,  ),  and  the  objective  f^ineticm  tekee  on  the  vmlue  22. 

We  have  here  deleted  all  n  reira  and  oolemne  (n  •  3)  opacified  eoeording  to 


Theorem  3 
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It  is  gen«rall^  posslbla  to  operate  with  n  fewer  rows  frora  the  start. 
We  say  in  fact  use  only  the  tableau 


-b 

I 

a'' 

». 

0 

w 

0 

In  order  to  apply  the  hybrid  dual  method  to  (7)t  bhe  algorltho  mtiat  be 
modified  PB  follows. 
ctsice  li  ^Aa  before.) 

Stare  A  If  Cj  <  0,  we  may  replace  Cj  with  Cj  -  [oj  I,  Cj  with 


CcJj,  a[j  with  ajj  -  [ajj],  and  ajj  with 


‘IJ 


for  1  •  We  do  not  conaider  ths  negative  of  the  colucm  in  which 

e 

c‘,  aopearr# 

If  c j  >0,  replace  Cj  with  Oj  ♦  c *  ,  Cj  with  Or 
with  ♦  a*j,  and  a*j  with  0  (for  1  - 

B.  £  Co  U  -  0.  the  row  of  tableau  (5)  cor.talnjjig  the  Jth  r<M  of 

V 

may  be  annexed  to  (7)t  and  the  method  of  Theorem  3  vaed  to  orlng 
the  Indicated  oolmn  anl  eome  row  into  the  required  fom  for  rcsoval. 

0<  The  eolution  is  obtained  in  the  same  fashion  as  for  the  regular 
hybrid  4lual  al^rltbm. 

We  naj  Justify  the  andifled  method  applied  to  tableau  (7)  by  observir.s 
that  it  is  newer  necessary  to  alter  the  'I^^n  ^  ***  sleet 

not  to  add  a  anltiple  of  any  of  the  other  rows  to  it«  The  adjustsents  of  the 
tableau  when  a  eompanent  of  is  negative  then  becoas  the  sane  ae  those 


•peolflad  by  Theor«a  1  In  creating  (S)n  The  adjustaente  wh«  a  eoaqponant  of 
is  poeltlve  slailarly  result  ty  adcnoel  edging  the  iaplieit  sodLstenoe  of 

-nXn 

Ae  tdth  TheoroB  1,  the  analysis  Is  stzlotly  correct  only  If  b  >  0, 

If  some  conpcnent  of  b  becooes  sero  it  nay  then  be  neeeseary  to  deviate 
from  the  preceding  to  maintain  lexicographic  negativity.  The  condltiona 
for  deviation  follow  the  ordinary  roles  for  applying  the  bound  esoalation 
method,  and  hence  do  not  require  additional  specification.  The  reeult  is 
that  one  or  more  of  the  c  rows  associated  with  the  ■etriz  nay  have 

to  be  ezpTdcltly  represented  in  ease  of  degenaraoy.  However,  we  suggest 
as  in  Section  111  that  the  departure  from  strict  leadeograpbio  negativity 
brought  about  by  ignoring  degonsrsoy  in  the  b  vector  relative  to  tho 
rows  of  will  not  nonaally  interfere  with  convergence.  ^ 

To  illustrate  the  preceding  ideas  more  thoroughly,  we  will  rework 
the  tuaaple  nroblma  roD.ved  earlier.  Using  ths  condenssd  repreeentatlon, 
ws  begin  with  the  first  tableau  for  Stags  2. 


(0) 

22/7 

4/7 

2/7 

6/7 

1 

0 

0 

0 

2 

-1 

-13/7 

3A 

5/7 

1/7 

0 

1 

0 

0 

0 

0 

5/7 

6/7 

3/7 

2/7 

0  . 

0 

1 

-1 

0 

_2- 

131/7 

1/7 

4/7 

5A 

0 

0 

0 

‘2 

-4 

Col.  3 

Col 

(i) 

-1 

-1 

0 

1 

0 

-3 

3 

2 

-1 

0 

1 

-8/7 

3/7 

2,-7 

-V7 

0 

1 

-1 

1 

0 

0 

2/7 

0 

-5/7 

6/7 

3/7 

2/7 

0 

0 

1 

.1 

0 

0 

3/7 

0 

U6/7 

-17/7 

-5A 

-1/7 

0* 

0 

-3 

1 

-1 

-4 

2/7 

-2 

9 


1.  If  laxicographio  nsgstivlty  is  lost  in  this  fashion  it  may  always  be 
regained  by  insertiug  a  ooluen  oonaiaUi^  of  all  nogative  eooffioisnte 
iBsediately  to  tho  right  of  the  b  vector,  with  a  aiffLelently  largo  nogative 
nuaoer  In  the  bottom  row  to  make  the  constraint  eompatiblo  with  any  optimal 
solution. 
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Vhe  two  colvs^'A  at  ths  end  of  the  tabljeau  indicate  the  colunna  to  replace 
colistn  3  co?.tznn  9  according  to  the  rule  of  the  oodifled  li^straction  Ac 
Snlectlog  the  now  eolvnn  3  as  the  next  colwi  to  put  in  the  d  fom, 

K«  obtaj,'!  the  following  tableau. 


(u)  -1 

i6 

o 

o 

1 

0 

-3 

3 

1 

*1 

*3/7 

V7  '1/7  •3/7 

0 

1 

2 

2 

0 

0 

.  ->/7. 

6A  3/7  2/7 

0 

0 

-1 

0 

0 

15V7 

23/7  -1/7  “3/7 

0 

0 

4 

2 

-<2 

•4 

(lil)  -1 

2 

0  0 

1 

0 

"3 

3 

1 

-1 

'5/7 

-9/7 

6A  ’3/7 

0 

1 

2 

2 

•  1 

0 

15/Y 

-1/13/7^ 

0 

1 

0 

22 

•2 

0  0 

0 

~1 

«2 

0 

2 

-4 

Col,  3  Col.  9 
0  1 

6A 

JdJL 

6/7  -3 


The  probl9tt  is  now  so?.v»<l,  and  the  final  solution  vs-^tor  is  obtained  by 
r.dding  (  2  0  0)  to  the  (0  -2  -4)  rector  on  the  botton  row,  yielding 
( 0  '  1  2-2  2  4}  as  baforS/ 

Two  points  of  contrast  may  te  noted  between  the  regular  hybiid  <iual 
algorithm  and  the  modified  method  relative  to  the  preceding  example  problem. 
First,  the  iixxlified  orithod  did  not  do  so  well  as  ttie  regular  nethod  in 
preparing  colwns  of  the  tableau  for  eleination  (though  the  siballer  dimcmdon 
of  (7)  aerred  adequate]^'  tu  c-omper.sate  for  this  in  the  present  Instance}, 

Secoial,  tits  modified  uethod  requited  an  additional  step  to  obtain  the  optdjb  l 
solutioQ.  It  may  be  obsetnred  that  with  the  regular  method  the  solution  wia 
obtained  >)  focusing  only  on  coliwn  4,  Following  this  dtratsgy  with  the  arxiifisd 
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Bsethod.  requires  jet  another  step  bsyronl  Uie  nuadber  needea  aboTSo  The  reaso. 
vhj  this  icight  bo  so  is  suggested  in  [3],  where  It  is  shown  that  It  Is  generailjr 
desirable  to  make  row  x  additions  with  the  bound  eeealation  Mthod  wheneeer 
pnaslbleo  Suppressing  the  aattix  in  effect  rules  out  the  chance 

for  aocM  of  these  additions.  Thus,  it  may  be  in  general  that  the  modified 
aathod  requires  more  steps  to  obt.'iln  an  optijsal  solution  than  the  regular 
methodi  perhr.pe  partialljr  offsetting  the  savings  derived  frcB  using  a 
uaaller  tableau.  Knreover,  the  initial  difference  in  tableau  else  may  be  a 
transltoxx  une  due  tc  the  superior  ability  of  the  regular  hybrid  dual 
algorithm  to  ellmirate  rows  and  colunns. 

ihere  is  another  wa>  of  working  with  a  substantially  smaller  tableau 
than  (5)»  while  still  retaining  the  ability  to  perfom  row  additioos 


relati  re  to  the  *'1^^  matrix.  The  key  la  to  consider  the  oolune  sssoeiat«<i 


with  c"  onq  at  a  time,  applying  the  nethod  of  Theoren  3  to  each  in  sueesssi  'n 
It  is  evident  triat  the  elensnt  d^  of  d  does  not  have  to  begin  tero,  but  maj, 
be  driven  to  sero  exactly  as  the  elsnents  of  Thus  each  oolum  (and  an  aisiciated 
row;  Is  eliminated  before  movli^  to  the  nett  (mises  the  aolutlan  is  obtaii  >d 
first).  In  thde  way  only  s  sin^e  row  associated  with  ttio  mstrix  net<;s 

to  be  iroludsd  In  the  tableau  at  any  stage,  bueh  a  prooedurs  corrsspoads 
to  the  steps  setusUy  taken  in  solving  the  eteaT>1s  problem  in  this  paper*  and 
seec«  to  be  about  mm  afructive  as  any  ether  method  for  solving  the  slmpls 
problems  on  which  we  have  tested  the  algorithm  to  date.  If  this  method 


were  used  consistently  a  fui^her  savinas  in  t^Me'ni  •ire  scoxd  bo  effected 
by  lottix^  the  initial  aatrix  be  ecual  to  the  A  natrlx  of  the  final 
linear  inogr-iimoing  tableau,  exceut  for  the  single  column  to  which  the  methed 


-21 


of  Theore.-i  3  was  to  be  applied^  olnce  none  of  tbe  operations  would 
affect  ccluotn  b  j  of  until  colunn  J  of  A  was  itself 

operated  on  directly  fay  Theorem  3»  not  all  of  the  eoluaia  of  rm;x(B-m> 
%«oula  need  to  be  given  explicit  representation  This  would  in  effect 
produce  a  constant  tableau  size  n  -  1  rows  and  n  -  1  oolumns  soaller  than 
(3)  until  the  last  coluaui  of  A^  vras  elijainated,  for  an  old  row  and  coluon 
would  be  dropped  each  tine  a  new  row  and  oolusm  required  explicit 
represenlstionc  A  possibility  that  iiii^t  be  worth  considering  would  be 
to  segregate  the  vector  tc  which  Theorem  3  was  to  be' applied  (along  with 
the  objective  function  vector until  it  was  reduced  to  the  indicated 
form  for  removal.  Tliereupon  the  transfomation  which  caused  the  reduction 
could  be  applied  to  update  the  reaainder  of  the  tableau. 

Though  we  have  not  exemplified  this  second  oethod  of  tableau 
reduction,  since  it  is  relatively  straightforward,  we  would  expect  it  to 
be  r^eferable  to  the  modified  metaod  based  on  table<'iu  (7>  in  m^ny 
cir>niiLStar.ces.  As  a  gmieral  rule,  however,  the  price  of  using  either  of 
the  saaller  tableaus  is  a  reduction  iJi  the  nuaber  of  alternatives  of 
choice,  >nd  we  have  yet  to  oiscovar  in  vdiich  cases  tMs  price  is  greater 
or  less. 

SooB  cuiBoer.ts  on  the  fera  of  the  A^  matrix  are  xn  order.  It  was 

shown  in  Ineorcx  1  that  this  uatiix  results  uloply  from  applying  the  bound 
SBcalation  method  to  ebcti  of  the  -nlunns  associated  with  the  A  matrix  of 

tableau  (6),  It  ie  interesting  to  note  that  the  resulting  constTulnts  are 
In  fact  "Goiaory  cuts"  of  the  type  used  in  his  first  integer  algorithm  By 
deriving  these  constraints  with  the  basic  transform t ions  prescribed  by 
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the  bound  escalation  method,  we  have  shovm  how  each  oiw  cay  not  only 
be  appended  to  the  tableau  as  In  (kmoxy’s  flret  algorithm  (in  addition 
to  the  constraint  from  which  it  was  derived),  but  abjr  In  fact  "replace" 
its  parent  constraint,  provided  our  rules  for  operating  on  ($)  are 
observed,*^  Moreover,  instead  of  using  the  linear  prograsaing  pivot 
process,  khieh  nay  require  the  generation  of  other  derivative  cuts,  we 
pro  reed  directly  with  the  bound  escaXatlcm  jsethod  in  conjunction  with 
Theorem  3.  shelving  off  poi  tioiis  cf  the  tableau  as  they  a  ttaln  the  qieeiiiel 
form 

Tc  a  rouji  Uoa  jT  the  compatational  perforaanec  of  Iho  I^-br^r.- 

duel  algorithm  we  have  solved  eXeven  eddlticnal  problwe  with  the  nethod 
For  eocparleon.  those  problems  were  also  solved  with  the  two  integer 
pro^rcmnalng  algorithms  developed  by  Ralfh  Gcmoiy*  The  ceope  of  the  teetin^ 
was  eotreoMly  limitedx  none  of  the  problems  was  larger  than  three  varlablrs 
and  three  ineqtialltiss.  No  reel  conelusione,  of  eouree,  can  be  drawn.  However, 
for  t’.  c  xTob?(rre  veauinwd,  the  hybrld<-duel  method  required  fewer  stepe  ii  laoet 
every  case  than  the  nmt'yT  of  pivots  required  by  the  all*>iiiteger  algorithm 

1.  For  eaditional  (non-Gonoiy )  cuts  wiiich  may  be  obtained  with  the 
transforaations  ^eeifled  in  tuid  which  alternately  may  be  used  as  a 
basis  for  the  initial  tableau  of  Stage  2  of  the  hybrid-dual  algorithm, 
see  C7J  An  interesting  feature  of  these  latter  cute  ie  thut  both 
ixMitive  and  negative  coefficients  may  ettaeh  to  the  varlablee  aaaoelated 
with  the  A  matrix  of  (o),  ae  contrasted  with  the  all  noonegative 
coefficients  of  A^, 


-5>- 


On  the  other  hand,  no  narked  difference  vae  obeerenbie  between  the  bgfbrld-^ 

e.\po»dth«  an..  oTdgltMLl  Integer  progreB  o-Ot:  Q«M'>odf  the  b^hrid 

dual  algorltha  doing  elightlj  better  on  eooe  of  the  problene  end  alightljr 

worse  on  others.  Both  nethods  soleed  8  of  the  11  probleao  in  froa  one  to 

# 

three  etepe  after  the  linear  prog<«aedng  solution  was  Outaiiiod,  urd  szlrzi 
romc^Lning  three  problene  within  six  steps  after  tlte  linear  prograunlng 
solut  An  edge  ahould  be  given  to  Uonory'e  first  integer  ;.rograaalng 

'-.Igor^^hu  fur  tlto  prc.-bldns  teatec,  hoMver,  because  of  a  analler  average 
tableau  size.^  Further  investigation  is  noeded  to  detemine  relative 
perforaances  for  prohlees  of  higher  diaenaion  and  eoeplexitj. 

It  might  be  noted  that»  because  of  the  sii^ilicity  of  the  problent. 
examined,  coiqjutatlonal  adv-'intagee  of  the  general  bound  escalation  aeU.od 
(in  contract  to  the  simpler  version  presented  in  this  pe^er)  were  not 
exploited  For  problacts  in  vdilch  tbs  D.  Q.  d  structure  may  be  manufacturec 
to  exhibit  certain  ohardcteristlce  (sec.  e.g,.  example  problv  2  of  [3i>« 
the  general  bound  earalation  method  by  itself  is  appreciably  superior  to  th«.' 
algcrithms  diecuaeed  here.  It  would  be  interesting  to  know  whether  its 
incorporation  into  the  hyhrld-<dual  algorithm  might  prove  sifflilarly  effective 
in  solving  etill  other  classes  of  probletoa  idiich  have  not  readlijr  yieldec 
to  either  of  Cioaory's  algorithms  or  to  the  bound  escalation  method  in  its 
independent  form. 


lo  The  second  procecure  outlinec  in  this  paper  for  reducing  the  dimension 
of  (5>  woulO  have  substantiall/  eliminated  thia  difference.  Ue  did  not 
apply  the  modlfsd  method  based  on  tableau  (7). 
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